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ABSTRACT 

[-] Mod e lling A modelling method having application to petroleum production, soil 
cleaning, etc. for optimizing faster and more realistically the displacement conditions, in 
a porous medium wettable by a first fluid (water for example), of a mixture of fluids 
including this wetting fluid, another, non-wetting fluid (oil for example) and a gas. 
£-l-The method comprises experimental determination of the variation curve of the 
capillary pressure in the pores as a function of the saturation in the liquid phases, 
modelling the pores of the porous medium by means of a distribution of capillaries with 
a fractal distribution by considering, in the case of a three-phase water (wetting fluid)- 
oil-gas m ixture for example, a stratification of the constituents in the pores, with the 
water in contact with the walls, the gas in the eanfre -center and the oil forming an 
intercalary layer, determination, from this capillary pressure curve, of the fractal 
dimension values corresponding to a series of given values of the saturation in the liquid 
phase, modelling the hysteresis effects that modify the mobile saturations of the fluids 
effectively displaced in the sample, that vary during drainage and imbibition cyclesj. 
modelling th e relative pormoabiliti e s dir e ctly in tho form of analytic expr e ssions 
depending on the various fractal dim e nsion values obtain e d and in accordance with the 
stratifi e d distribution of tho various fluids in th e poroo, and using a porous medium 
s imulator in order to dotormino, from th e r e lative permeabilities, tho optimum 
displacem e nt conditions of tho fluids in th e porous m e dium. 
■ Application : p e troleum production, soil d e poUution, etc. 
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Field of the Invention 



FIELD OF THE D^JVENTION 



rooon 



The present invention relates to a method for modelling two-phase or 



three-phase flows in a porous medium, in drainage and imbibition. It is based on a 
fractal representation of the porous medium and on an original approach for handling 
phenomena linked with hysteresis (change in the direction of variation of the 
saturations). 

Th e m e thod according to the invention is applicable in many fields whoro fluid 
flows in porous m e dia a r e t o b e m od e ll e d i n o rd e r too ptimiz e t h e o onditions o f t he 
displacem e nt th e r e of in drainag e and imbibition. Examples of fields of application ar e : 

a) d e velopm e nt of an oil r e servoir and notably enhanc e d hydrocarbon production by 
injection of fluids, using for e xampl e alt e mat e inj e ctions of liquid and gas slugs (a 
method referr e d to as WAG). It constitut e s an advantageous tool allowing res e rvoir 
engineers to study also well productivity and injectivity probl e ms ; 

b) soil dopoUution and notably d e poUution of industrial sites by inj e ction of 
substanc e s such as surfactants in pollut e d layers; 

c) cleaning of filt e r cak e s by displacement of th e substanc e s r e tain e d th e r e in; 

d) wood drying ; 

o) optimization of chemical r e actions for e xampl e by displac e m e nt of r e action 
products in a catalyst mass in ord e r to incr e as e th e surfaces of contact, e tc. 
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Tho method according to th e invention is directly applicabl e by r e s e rvoir e ngin ee rs 
in ord e r to dotormino, for oxamplo, tho most suitable enhanc e d r e covery m e thod to b e 
appli e d to an xmdorground hydrocarbon reservoir. Tho method con also servo within tho 
scop e of industrial sit e depoUution op e rations for example, 

BACKGROUND OF THE IhJVENTION 

Description of the Prior Art 

1) Experimental studies 

r00021 ^Experimental determination of the relative permeabilities of a porous 

medium wherein a multiphase fluid flows is not easy. Measiuing operations are usually 
simplified by considering that one of the phases is inraiobile in a state of irreducible 
saturation. 

r00031 ^The values are for example acquired by means of a well-known 

experimental method referred to as "steady state" for determining relative 
permeabilities, which consists in allowing allows a three-phase fluid to flow with 
imposed flow rates between the phases. The relative permeabilities expressed as a 
function of the two saturations are calculated by applying Darcy's law to each phase. It is 
not an established fact that the relative permeability measurements obtained by means of 
this method are really representative of the fluid displacements and, in any case, they 
take a long time because, at each regime change, one has to wait for a state of 
equilibriimi. 

r00041 Another known method consists in carrying carries o ut laboratory tests in 

order to determine measurement tables (as shown in Fig.l) relating the relative 
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permeabilities and the saturations for each pair of fluids of the three-phase mixture. By 
adjusting experimental production curves, one tries to progressively adjust the three- 
phase relative permeabilities. These data tables are then entered into an Athos® type 

numerical simulator which computes the fluid productions. This method bem^is based 
on the prior acquisition of many experimental measurements progressively adjusted by 
calibrationr4 t and takes a long time. 

2) Relative permeability models 

rOOOSI The known empirical model referred to as Stone's model allows, by 

empirical correlations, to predict data relative to a three-phase flow from data 
corresponding to a two-phase flow. It is valid only in case of a high water wettability 
and it is generally considered to b e poorly pr e dictiv e a poor predictor . 

fOOOSI ^There are two known types of physical models for modelling three-phase 

flows, based on capillary pressure curves. The capillary pressure curves are connected 
with a saturation (for example that of the mercury injected) and a pore radius from 
which the mercury stops for a given injection pressure, determined by Laplace's law, 

Pinj.= ^. 
r 

r00071 A first porous media representation model is described by : 

- Burdine, N.T. : "Relativ e perm e ability calculations from por e siz e distribution 
data ^" Relative Permeability Calculations from Pore Size Distribution Data". Trans 
AIME (1953), Vol.198, or by 
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- Corey, A.T. : "The int e rr e lation b e tw ee n oil and gas r e lativ e p e rmeabiliti e s"] [The 
Interrelation between Oil and Gas Relative Permeabilities", Prod. Monthly (1954), 
Vol.19, 38. 

rOOOSI According to this model, the porous medium is represented by a bundle 

of cylindrical capillaries with a radius distribution given by the capillary pressure curve 
obtained by mercury injection. The permeabilities are obtained by applying Poiseuille's 
law to the flow of fluids in these capillaries. 

r00091 ^This model is based on the representation of the porous medium as an 

assembly of capillaries with different radii. The relation between the volume and the 
radius of the pores is given by the value of the slope of the pseudo-plateau. The three 
fluids are supposed to share the capillaries between them, the wetting fluid (water) 
occupying the smallest ones, the least wetting fluid (gas) the largest ones and the third 
fluid (oil) a zone with intermediate-size pores. It is not possible to describe the 
interactions between the fluids because, in such a model, they flow through separate 
channels. Finally, this model can be useful only if the pseudo-plateau covers a wide 
range o f s aturations. According to this model, the three phases of a three-phase flow 
move in different capillaries and there is no interaction between them. 

[00101 Another known physical porous medium representation model is 

described 
by: 
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- de Gennes, P.G. : "Partial Filling of a Fractal Structure by a Wetting Fluid", 
Physics of Disordered Materials 227-241, New York Plenum Pub. Corp. (1985), taken 
up by 

- Lenormand, R. : "Gravity Assisted Inert Gas Injection : Micromodel Experiments 
and Model based on Fractal Roughness", The European Oil and Gas Conference 
Altavilla Milica, Palermo, Sicily (1990). 

rOOin According to this model, on e considers that the inner surface of the pores 

is considered isotropic and has a fractal character, and it can be modelled as a "bunch" 
of parallel capillary grooves so that the pores exhibit a fractal cross section. The cross 
section of each pore is constructed according to an iterative process (Fig.l). The half- 
perimeter of a circle of radius Rq is divided into r| parts and each of these r\ parts is 
replaced by a semi-circle or groove. At each stage k of the process, new semi- 
circular grooves of radius R^ and of total section A^ are created. 

r00121 ^The fractal dimension DL o f t he c ross s ection attheendofs tage k i s 

related to the number of objects Nk generated with the given scale 1^ by the relation : 

r00131 ^The fractal dimension can be deduced from a mercury capillary pressure 

curve according to the following procedure. Mercury is injected into a porous medivun 
with an injection pressure that increases in stages. Laplace's law allows to d e duc e 
deducing the pore volume, knowing the volume of mercury injected for a given capillary 
pressure and on e can construct the drainage capillary pressure curve relating the 
injection pressure to the amoimt of mercury injected and the curve relating the 



6 



proportion of the total volume occupied by the pores and the size of the pores can be 
constructed . In cases where a wetting liquid is drained from the porous medivmi such as 
water by gas injection, the correlation between the gas- water capillary pressure and the 
saturation of the wetting phase is given by : 

1 



r0014l ^The experimental results readily show that the values of the gas-water 

relative permeabilities expressed as a function of the three saturations, obtained from the 
expressions given by the known models and the phase distribution modes in the 
structure of the pores, are far from the measured values and therefore that the models 
concemed prove to be too simplistic to represent the complex interactions that take 
place between the fluid phases. 

r001 51 P atent FR French Patent 2 ,772,483 (US-6,02 1,662) describes a modelling 
method for optimizing faster and more realistically the flow conditions, in a porous 
mediiun wettable by a first fluid (water for example), of a mixture of fluids including 
this wetting fluid and at least another fluid (oil and possibly gas). This method involves 
modelling the pores of the porous medium by a distribution of capillaries with a fractal 
distribution considering, in the case of a three-phase water (wetting fluid)-oil-gas 
mixture for example, a stratification of the constituents in the pores, with the water in 
contact with the walls, the gas in the c e ntr e center and the oil forming an intercalary 
layer. It comprises experimental determination of the variation curve of the capillary 
pressure in the pores as a function of the saturation in the liquid phases, from which the 
fractal dimension values corresponding to a series of given values of the satiu-ation in 
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the liquid phase are deduced. It also comprises modelling the relative permeabilities 
directly in_a form of analytic expressions depending on the various fractal dimension 
values obtained and in accordance with the stratified distribution of the different fluids 
in the pores. A porous medium simulator is used from these relative permeabilities to 
determine the optimum conditions of displacement of the fluids in the porous medium. 

r0016l ^The hysteresis phenomenon relates to the variations in the petrophysical 

properties (relative permeabiUties, capillary pressure, resistivity index, etc.) observed 
according to whether measurements are performed during drainage or imbibition (these 
modes respectively correspond to a saturation increase and decrease of the non-wetting 
phase). This phenomenon must therefore be taken into account to obtain representative 
relative permeability values. 

r00171 ^The prior art concerning hysteresis effects in two-phase and three-phase 

media is described for example in the following publications : 

- Land C.S. : « Calculation of imbibition r e lativ e p e rm e ability for two and thr e e 
phas e flow from rock properti e s » "Calculation of Imbibition Relative Permeability for 
Two and Three-Phase Flow from Rock Properties", Trans AIME 1968, Vol.243, 149, 

Larsen J. A., Skauge A. : « Methodology for numerical simulation with oyolo 
d e pend e nt relativ e p e rm e ability », "Methodology for Numerical Simulation with Cvcle- 
Dependent Relative Permeability", SPEJ. June 1998, and 

Carlson P.M. : « Simulation of relative permeability hysteresis to the non 
w e tting phas e )) "Simulation of Relative Permeability Hysteresis to the Non- Wetting 
Phase", SPE 10157, ATCE, San Antonio, Texas, 4-7 Oct. 1981. 
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[00181 Figur e Fig. 6 t ypically s hows the course of the two-phase permeability 

curves Knv resulting from drainage up to irreducible saturation in wetting fluid (M), then 
imbibition up to residual saturation in non-wetting fluid (MM). The hysteresis 
phenomenon occurs at two levels. At equal saturations Sg, difTerent numerical values are 
obtained and the end point reached is an unknown that depends on the cusp point SgM 
from which the displacement mode is changed. This phenomenon is usually attributed to 
the trapped non- wetting fluid fraction. At equal saturations, w e thus do not hav e the 
same quantity of mobile flui d is not obtained, which distorts the flow characteristics. 

r00191 ^Practically all the models taking account of hysteresis effects involve 

Land's semi-empirical relation : — — = (1) 

where Cl represents Land's constant. This relation relates the initial saturation Sgi to the 
residual saturation Sgr in non-wetting fluid, in order to evaluate the saturations in 
trapped and free non-wetting fluid. Assuming that this relation is valid whatever the 
saturation, it is applied to determine the intermediate mobile fractions during 
displacement. In the case of two-phase flows, associating this relation with permeability 
models provides satisfactory results. 

r00201 Li the case of three-phase flows, the hysteresis of the relative 

permeabilities Krg takes on a particular form. A displacement hysteresis is 
experimentally observed, as it is the case with two-phase flows (extent of the direction 
of variation of the saturations), as well as a cycle hysteresis since the permeabilities 
depend on the saturation record. In Figur e F ig. 7, the relative permeability curves Krg 
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corresponding to a first drainage and imbibition cycle (Dl and II respectively) are 
distinct fi"om the corresponding curves (D2, 12) of a second cycle. 

r002n The relative permeability model developped by Larsen et al takes these 

two forms of hysteresis into account. Starting from an approach combining Stone's 
model in parallel with Land's formula and Carlson's interpolation method, an approach 
where only the displacement hysteresis is taken into accoxmt, Larsen et al have 
introduced an empirical reduction factor that is a function of the water saturation, which 
allows to approximate to the permeability reduction of the gas associated with the cycle 
hysteresis. 

SUMMARY OF THE INVENTION 

r00221 The modelling method according to the invention allows to optimiz e 

festef -faster optimizing and more realistically the displacement conditions, in a porous 
mediimi wettable by a first fluid, of two-phase or three-phase mixtures including this-the 
first wetting fluid and at least a second, non-wetting fluid. It therefore provides 
operators with a more reliable tool for evaluating notably the best displacement modes 
of the fluids in the porous medium, in drainage and imbibition. It is based on a fractal 
representation of the porous medium with modelling of the pores by a distribution of 
capillaries with a fractal section, considering a stratified distribution of the fluids in the 
pores, the wetting fluid being distributed in contact with the walls and aroxmd the 
second fluid (or around the other two in case of a three-phase mixture). 
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f00231 The method according to the invention is applicable in many fields where 

fluid flows in porous media are to be modelled in order to optimize the conditions of the 
displacement thereof in drainage and imbibition. Examples of fields of application are : 

a) development of an oil reservoir and notably enhanced hydrocarbon production by 
injection of fluids, using for example alternate injections of liquid and gas slugs (a 
method referred to as WAGV It constitutes an advantageous tool allowing reservoir 
engineers to study also well productivity and injectivity problems : 

b) soil depoUution and notably depoUution of industrial sites by injection of 
substances such as surfactants in polluted layers: 

c) cleaning of filter cakes by displacement of the substances retained therein: 

d) wood drying : 

e) optimization of chemical reactions for example by displacement of reaction 
products in a catalyst mass in order to increase the surfaces of contact, etc. 

[0024] The method according to the invention is directly applicable by reservoir 

engineers in order to determine, for example, the most suitable enhanced recovery 
method to be applied to an underground hydrocarbon reservoir. The method can also 
serve within the scope of industrial site depoUution operations for example. 

[00251 ^The method comprises in combination : 

- experimental determination of the variation curve of the capillary pressure (Pc) in 
the pores of a sample of this porous medium in the presence of a wetting fluid and of at 
least one non- wetting fluid (by injection of mercury in a sample placed under vacuum 
for example); 
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determination, from this capillary pressure curve, of the fractal dimension values 
corresponding to a series of given values of the saturation in liquids; 

modelling the hysteresis effects that modify the mobile saturations of the fluids 
effectively displaced as a function of the number of drainage and imbibition cycles 
undergone by the sample, involving different non-wetting fluid trapping or untrapping 
constants according to whether a drainage stage or an imbibition stage is carried out ; 

- modelling the relative permeabilities directly in the form of analytic expressions 
depending on the various fractal dimension values obtainedj; and 

- entering the relative permeabilities in a porous medium simulator and determining, 
by means of this simulator, the optimum displacement conditions of the fluids of the 
mixture in the porous medium. 

[00261 ^The method is applied for example for determining displacements of 

fluid mixtures comprising a first wetting fluid, a second, non-wetting fluid and a gas, 
considering a stratified distribution of the fluids in the pores, the wetting fluid spreading 
out in contact with the walls, the gas occupying the c e ntr e center of the pores and the 
second fluid being distributed in the form of an annular layer in contact with both the 
gas and the first fluid. 

f00271 ^The method can notably be applied for determining, by means of a 

reservoir simulator, the optimum characteristics of substances added to wetting fluid 
slugs injected in a formation alternately with gas slugs, in order to displace 
hydrocarbons in place, or those of a fluid injected into the soil in order to displace 
polluting substances. 
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f00281 Modelling of phenomena by means of the present method has many 

advantages. It allows a better correspondence with the results obtained in the laboratory 
because the physical phenomena are better taken into account. The results of the model 
are therefore better in case of a scale change for example, for modelling an application 
in an operation field. 

[00291 The calculating time is reduced by comparison with the time required 

when tables are used as in the prior methods. Fractal type modelling can better deal with 
the hysteresis effects encovmtered when using WAG type injection processes. 

[00301 The results of the method can furthermore be perfectly integrated into 

many reservoir simulators : 3D, heterogeneous, composition simulators, etc. 

[00311 Exploitation of the results by application softwares is facilitated. It is no 

longer necessary to perform risky interpolations as it is the rule when working from 
discrete values of the result tables in order to draw isoperms for example. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[00321 ^Other features and advantages of the method according to the invention 

will be clear from reading the description hereafter of a non limitative example, with 
reference to the accompanying drawings wherein : 

[00331 [-] Fi|gur e F ig. 1 illustrates, in the form of a table, the connections 

existing for a three-phase mixture between the experimentally obtained relative 
permeability values of a fluid and the saturations for two of the three fluids^; 

[00341 [-] Figure Fig. 2 is a fractal representation of a pore^; 
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f00351 [-1 Figur e F ig. 3 diagrammatically shows the distribution of the phases of 

a three-phase fluid in a fractal pore with the wetting fluid W in contact with the wall, the 
gas phase G spread over the greatest part of the volume of the pore (radius Rq)^ the oil 
O being a layer between the wetting fluid and the gas^; 

f00361 [-1 Fijguro Fig. 4 shows, as a function of the water saturation, the capillary 

pressure c urve o f a sample of Vosges sandstone for example, whose local slope S is 
used to determine the fractal dimension of the poresj; 

r00371 [-] Figur e Fig. 5 shows, as a function of time, the production curves 

obtained experimentally for gas (G), for water (W) and for oil (O), compared with the 
equivalent curves obtained by simulation by means of the method according to the 
invention^; 

f00381 ['} Figur e s F igs. 6a, 6B show the effects of the displacement hysteresis 

affecting the relative permeabilities of the wetting fluid and of the gas, and K^g 
respectively^^; 

[00391 [-] Figur e Fig. 7 shows the combined effects of the displacement and 

cycle hysteresis observed experimentally on the gas phgise^; 

r00401 F- l Figure F ig. 8 shows the key curves used for modelling the untrapping 

constant?; 

[00411 [-] Figure Fig, 9 shows the result obtained with the model within the 

scope of a WAG type injection?; and 

[00421 r - 1 Figure Fig. 10 shows the validation of the model on a WAG type 
experiment for the recovery curves of the three phases. 
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DETAILED DESCRIPTION OF THE METHOD INVENTION 

r00431 ^The method according to the invention allows to d e t e rmi fte -determination 

of the three-phase relative permeabilities of porous media by using a fractal type model 
of the porous medium, on the basis of an approach described by : 

- Kalaydjian, FJ-M et al. : "Three Phase Flow in Water-wet Porous Media : 
Determination of Gas-oil Relative Permeabilities under Various Spreading conditions", 
68th Ann. Tech. Conf. and Exh. of the SPE, Houston, Texas, 1993. 

[00441 ^The method according to the invention comprises, as mentioned above, 

modelling the flow of the phases with distribution of the phases in the fractal structure 
of each pore. In the case of water and oil flowing in a water wet porous structure, the oil 
flows into the volmne of the pore. In the case of a three-phase flow, there is a stratified 
distribution, the water, which is the wetting fluid, flows along the walls of the pores, the 
gas circulates in the volume of the pore and the oil flows between the gas and the water. 
The saturations are calculated as the relative surface area in a cross section occupied by 
each of the fluids. 

r00451 ^At equilibrium, all the grooves having a radius greater than Rk, which is 

given by Laplace's law Pc = 2 y/Rk, are occupied by the gas, and the smallest tubes by 
the two other fluids (water and oil). The wetting fluid saturation is thus expressed as the 
fraction of the area of the occupied tubes. 

[00461 Calculation of the fraction of the area of the capillaries occupied by the 

water for all the radii between Rj^ and Roo leads to the following expression : 
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and, as Pc = 2 y/Rk, the correlation between the capillary pressure and the saturation of 
the wetting phase is given by : 



where Sw is the saturation of the wetting phase. 

[00471 The graphical representation of this correlation in a log-log diagram is a 

straight line starting from the point (Sw, Pc) corresponding to the largest capillary of the 
fractal structure with a radius Rq. One may suppose that : 

- radius Rq, which is first invaded when mercury is injected (Fig.4), corresponds to a 
saturation value of the order of 1/r = lO^. Each segment of the capillary pressure curve 
is a part of a line starting from Rq (assumed to be the same for all the different 
segments), corresponding to the aforementioned correlation Pc, Sw. Each line has a given 
slope, a fractal linear dimension can be associated therewith. The values of the slope 
range from 
-1.5 to -3.3 as shown in Fig.4, which leads to values of the fractal linear dimension Dl 
between 1 .3 and 1 .7; 

- each domain is reached by the mercury for saturations corresponding to the place 
where Rq is found on each line. 

f00481 ^The saturation of the two liquids when the gas is present in the pore is 

calculated as explained above for a phase : 
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assuming that the two hquids occupy the tubes whose radii are less than or equal to 
and the gas, the centre of each pore. The oil saturation is the relative area of the cross 
section of the grooves occupied by oil, whose radius is less than or equal to K\^. 

Hvsteresis modelling 



phase, of the fraction circulating therein and consequently systematic estimation of the 
saturations corresponding to the stagnant fractions. This must be done for the two cases 
studied, for example drainage of water and oil by gas, and water imbibition. 



talc e s account of accounts for t he hysteresis directly at its origin, irer-thaUs^at the level of 
the non-wetting phase trapping and imtrapping phenomena. In Ficau- e Fig. 8. if the 
curves could be deduced exactly by translation, this would mean that part of the gas has 
been trapped during the secondary drainage and does not take part in the flow. In fact, 
this is not the case but the fact that, at equal gas saturation, the permeability is lower 
during D2 than during II means that the trapped gas fraction is larger during drainage. In 
other words, the non-reversibility of the permeability curves can be explained by a non- 
reversibility (hysteresis) between the trapping phenomenon and the imtrapping 
phenomenon. 



r00491 



Calculation of the relative permeabilities requires determination, for each 



roosol 



The original feature of the modeUing procedure propos e d h e r e is that it 



r0051l ^Land's formula (Equation 1) is thus kept, but an untrapping constant, 

valid during the drainage stages, different from Land's constant valid during the 



17 



imbibition stages, is introduced. It thus all comes down to modelling the evolution of 
the imtrapping constant during the cycles. 

[00521 ^The characteristics to be taken into account in the formulation are as 

follows. We consid e r It is considered t hat Cp = Cl during all the trapping stages because 
wo oro o f being directly in Land's conditions of application where a single constant is 
enough to describe the phenomenon. When is low at the end of the drainage process, 
Cd is close to Cp. This simply means here that the gas is easier to imtrap when a small 
amount thereof is present in the sample, because eHe-it_can th e consider b e considered 
that it is present in large pores that can be readily reconnected during re-injection. When 
Sgr increases, Cd becomes greater than Cl, which shows a less efficient untrapping 
process, 

[00531 ^For high gas saturations, a certain reversibility can be encountered, but a 

trapped gas fraction remains inacessible in small pores, A low single curve Krg is thus 
reached, which corresponds to a pseudo two-phase case where the oil is no longer 
mobile, 

f00541 ^According to this representation, based on experimental observations and 

working hypotheses, the value of Cd passes through a maximum since the untrapping 
constant is equal to Cl when Sg^ is low and when it is maximum, when the low-mobility 
curve is described. The following expression allows to r e conoilo reconciUation of all the 
previous aspects : 



— K 

rg rg min 



K 



D\ 
rg 



(2) 



'gr2 J 
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roossi 



,This formulation comprises several parameters : 



Cdm ' it takes into accoimt the difference between trapping and imtrapping, 
Cl :Land's constant, 
K,Bmin • low-mobility curve, 
E : calibration parameter, 

K'rg: value of the relative permeability to gas at the beginning of the previous imbibition, 
KP\g : value of the relative permeability to gas on the first drainage curve for the gas 
saturation corresponding to K^g. 

r00561 ^Whatever the order and the nature of the cycle considered, relation (2) 

and Land's relation (1) allow to d e tenm ne- determination of the trapped and mobile 
saturations by means of Land's formula while taking into account the three-phase 
character of the hysteresis (displacement and cycles) without using empirical reduction 
factors. 



5^ = - 5 



g 




^iS^-S^r+^iS^-S^,)] (3) 



and Sg = Sg, + Sgf 



(4) 



C is equal to Cl or Co according to the displacement mode. 



Sgf : free gas saturation. 



Sgi : trapped gas saturation. 
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[0057] Calculation of the relative permeabilities 

A) Relative permeabilities of liquids 

Application of Poiseuille's law to each capillary of the bundle for the phase that 
occupies it -the capillary allows to calculate the water and oil relative permeabilities 
(Fig.5). 

r00581 Other experimental studies (Larsen et al) have shown that there is a 

relation allowing to relate the residual oil saturation during the cycles to the trapped gas 
fraction. 

Sor ~ (Sor)sgt-0 *~2^Sgi 

(Sor)sgt=o represents the saturation in residual oil left in place in the medium before gas is 
trapped. 

r00591 ^If w e consid e r i t is considered o nly the circulating fraction which 

contributes to the hydraulic conductivity, the relative permeabiHties for water and oil are 
expressed as follows : 

Ko = K^i2Ph)[(S, + S^y - (S^ + (1 - a)S^ + (S^^h^^y] 
K^ = (S^-R(l-a)S^y-SC, 

R is the reduction factor linked with the trapping of the non-wetting phase. 

f00601 ^In these expressions, it is useful to mention that : 

- the irreducible water saturation Swi is assumed to be stable. 
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- the size range of capillaries occupied by the mobile oil is calculated as the difference 
between the sizes of the capillaries occupied by the two liquids with the total liquid 
saturation Sl = Sq+S^ and those of the capillaries saturated in water and stagnant oil, 

- Kro(2ph.) is the value of the oil relative permeability determined by an imbibition test 
with water and oil. When only the water and oil phases are present and since the sample 
tested is water wet, the oil will flow through the section of the pore exactly as the gas in 
a three-phase flow. 

B) Gas relative permeability 

rOOSn Since gas is a non-wetting phase, it occupies the central space of the pore 

and it spreads towards t he p eriphery t hereof a s t he g as s aturation i ncreases, h owever 
without coming into contact with the solid wall (Fig.6). One consid e rs I t is considered 
that the gas circulates in a single pore whose radius Rg is given by the relation : 

Rg = R0+R1+R2+ ... +Rk, 
the gas permeability being then given by : 

Kr,=K^^i\-is,+s^,ry (5) 

with a = — ^ — , Dl being the linear fractal dimension of the porous medium and Sl 
2-Dl 

the total liquid saturation equal to 1-Sg. 

f00621 ^The relative permeability model allowing to calculat e calculating K p^y, 

and K^g is installed in a simulator such as ATHOS® or GENESYS©. This allows 
to calibrat e calibration of experiments carried out in the laboratory and to optimiz e 
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optimizing t he conditions to be satisfied in order to displace petroleum fluids in place in 
a reservoir, either by gas injection or by alternate injection of water and gas slugs (a 
method referred to as WAG), by taking account of the pressxire and temperature 
conditions prevailing at the production depth. 

Validation 

[00631 The method according to the invention has been validated by means of 

various types of experiments : 

gas was injected into porous media containing water and oil, under various 
conditions. It can be seen in Fig.5 for example that a very good accordance is obtained 
between the production curves of the three phases (water, oil, gas) obtained 
experimentally and those predicted by the reservoir simulator fed with the data obtained 
in accordance with the method, 

- gas and water were also alternately injected (WAG type injection). Fig. 10 shows 
that an excellent accordance is also obtained in this case for the three phases throughout 
the experiment. 



